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ABSTRACT 

We present a new model for the coronal structure of rapidly rotating solar-type 
stars. The presence of prominences trapped in co-rotation 2 to 5 stellar radii above 
the stellar surface has been taken as evidence that the coronae of these stars must 
be very extended. The observed surface magnetic fields, however, cannot contain X- 
ray emitting gas out to these distances. We present an alternative model: that these 
prominences are trapped in long thin loops embedded not in the X-ray corona, but 
in the wind. Above coronal helmet streamers, oppositely-directed wind-bearing field 
lines reconnect to form closed loops which then fill up with gas that was originally 
part of the wind. We demonstrate that static equilibria exist for these loops at a 
range of pressures and temperatures. The maximum loop height falls as the rotation 
rate increases, but rises as the loop temperature decreases. For a solar-mass star with 
rotation period 0.5 days, whose X-ray corona extends IR^ above the surface, loops at 
temperatures of lO^'K can extend out to 5K^,. 
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1 INTRODUCTION 

Since the days of the first X-ray satellites, the X-ray emission 
from solar-like stars has been a puzzle. For main-sequence 
stars it is a clear function of rotation rate, rising steeply 
for stars rotating faster than the Sun until rotation rates of 
about 15 — 20kms~^ where it reaches a plateau at Lx/Ltoi = 
10~^ IVilhu 198411 . Beyond this rotation rate is the "satu- 
rated" regime where the X-ray luminosity is independent of 
rotation rate. This "saturated" behaviour persists until ro- 
tation rates of about vsini > lOOkms"^, where the X-ray 
luminosity begins to decrease again. This regime is referred 
to as "supersaturated" JProsser et al. 19961 IRandich 1998ll . 

The initial increase seems linked to the enhanced dy- 
namo activity that accompanies increased rotation, but the 
causes of the so-called saturation and supersaturation are 
still a matter of debate. While dynamo activity within the 
convective region may saturate when the back-reaction of 
the fluid becomes significant, the X-ray emission produced 
when this field escapes into the corona is only an indirect 
measure of this process. Levels of X-ray emission may be 
determined by the physics of the corona as much as by 
the physics of sub-surface dynamo activity. Saturation of 
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X-ray emission may represent the inability of the surface, or 
the corona, to accommodate any more flux ( Vilhu 1984J , an 
inability of the energy release process to extract more en- 
ergy from the magnetic field, or centrifugal stripping of the 
corona at high rotation rates llJardine fc Unruh 199911 . 

A further puzzle is the nature of this enhanced emission 
from the rapid rotators. While early observations showing 
little rotational modulation of the X-ray emission indicated 
that the emission (and hence the coronae) must be very 
extended, some observations suggested that the emission 
came from compact (solar-like) regions at high latitudes 
on the stars where they did not sufl'cr from rotational self- 
eclipse l|S ingh, White fc D rake 1996 Siarkowski ct al. 19%| 
|Giampc^^^^iI^^96| ~|^ffrics 1998J. In this case, how- 
ever, a high density was required to explain the high 
X-ray luminosity. In the last few years, however, spec- 
tra obtained with FUSE, Chandra and XMM-Newton 
have made a significant impact on this problem by 
revealing high coronal de nsities (|Dupree et al. 1993 
Schrijver et al. 1995| [Brickhouse fc Dupree 1998 



Audard et al. 20011 IMewe et al. 20011 [ Young et al. 2001 
'Gfidel et al. 2001" "Sanz-Forcada, Maggio fc Micela 2003 



iSanz-Forcada, Brickhouse fc Duprcc 2003 ) . A critical Bep- 
poSAX observation of two flares on the rapid rotator AB 
Dor showed no rotational modulation of the flare decay 
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phases, which lasted for more than one stellar rotation 
l |Maggio et al. 2000| |. Flare modelling suggested the flaring 
regions were compact (0.3i?*) and therefore must have 
been located at high latitudes. The ability to locate emit- 
ting features is improving rapidly with the beginnings of 
"X-ray Doppler imaging" and further eclipse observations 
HBrickhouse, Dupree fc Young 2001| IGiidel et al. 200311 . 
Most recently, simultaneous X-ray spectra and sur- 
face magnetograms (obtained through Zeeman-Doppler 
imaging) have made it possible to couple the surface 
field to the coronal emission using field extrapolation 
techniques ([Jardine, Collier Cameron fc Donati 2002| 

[Hussain, Jardine & Colher'Cameroir'20'Ol' 
Hussain ot al. 2002 Jardine ct al. 2002 I . A consistent 
picture appears to be emerging of the coronae of these rapid 
rotators. They are compact and highly structured, with 
discrete emitting regions often situated at the high latitudes 
where Doppler imaging shows many of the large star spots 
to be located. The coronal densities are typically higher 
than on the Sun, giving the very high X-ray luminosities. 

One type of observation, however, appears to be 
quite inconsistent with this picture. On almost all these 
rapid rotators with He, observations, slingshot promi- 
nences are detected ([ Collier , Cameron fc Robinson 1989al 
ICoUier Cameron fc Rob inson 19893 

ICoUier Cameron fc Woods 19921 IJeffries 19931 

[ Byrne, Elbe fc RoUeston 1996[ lEibe 19981 

barnes et al. 20001 IDonati et al. 200011 . These are ob- 
served as transient absorption features that move through 
the Ha line profiles. In many instances these features 
re-appear on subsequent stellar rotations, often with some 
change in the time taken to travel through the line profile. 
These features are interpreted as arising from the presence 
of clouds of cool, dense gas co-rotating with the star and 
confined within its outer atmosphere. As many as six may 
be present in the observable hemisphere. What is most 
surprising about them is their location, which is inferred 
from the time taken for the absorption features to travel 
through the line profile. Values of several stellar radii from 
the stellar rotation axis are typically found, suggesting 
that the confinement of these clouds is enforced out to 
very large distances. Indeed the preferred location of these 
prominences appears to be at or beyond the equatorial 
stellar co-rotation radius, where the inward pull of gravity 
is exactly balanced by the outward pull of centrifugal 
forces. Beyond this point, the effective gravity (including 
the centrifugal acceleration) points outwards and the 
presence of a restraining force, such as the tension in a 
closed magnetic loop, is required to hold the prominence 
in place against centrifugal ejection. The presence of these 
prominences therefore immediately requires that the star 
have many closed loop systems that extend out for many 
stellar radii. 

While these observations may lend support to the hy- 
pothesis that these stars have coronae that are very ex- 
tended, it does not imply that they are smooth. Of all the 
multipolar components of a coronal field, the dipole compo- 
nent falls of most slowly with height, and so it is natural to 
suppose that a very extended corona will resemble a dipole 
field. A smooth and symmetric field like this, however, would 
give prominence locations which were, by symmetry, all in 
the equatorial plane and uniformly spread around the star 




Figure 1. A schematic diagram of the magnetic field in the equa- 
torial plane of the star as seen by an observer looking along the 
rotation axis. Wind-bearing, open field lines emerge from between 
regions of closed field that are bright in X-rays. Above these closed 
field regions, helmet streamers form. It is above the cusps of these 
helmet streamers that prominence-bearing loops may form. 

ll Jardine et al. 20011 . Prominences in the equatorial plane of 
a star are unlikely to be observed, however, unless the incli- 
nation of the stellar rotation axis is close to 90° . For a star 
such as AB Dor, whose rotation axis is inclined at some 60° 
to the observer, prominences that are located at around 3i?* 
from the rotation a^cis must be located at high latitudes in 
order to transit the disc. The large number of prominences 
observed (typically 6 in the observable hemisphere for AB 
Dor) suggest that there is significant azimuthal structure in 
the field supporting the prominences even at several stellar 
radii from the rotation axis. 

This high degree of complexity in the coronal 
field is supported by the degree of complexity seem 
in the surface magnetograms acquired by Zeeman- 
Doppler imaging and by traditional Doppler imaging 
IIDonati fc Collier Cameron 19971 IDonati et al. 199911 . Mul- 
tipolar flux is seen at all latitudes on AB Dor, even at very 
high latitudes. One region where it is not possible to deter- 
mine the field polarity is, however, very close to the rotation 
pole where the Zeeman signature of often suppressed by the 
presence of the polar spot that is typically found on rapid 
rotators ( Strassmcicr 1996J . By modelling the effect on the 
coronal structure of placing a unipolar field in this dark polar 
cap, jJardine et al. (2001)| jMcIvor et al. (2003)| have shown 
that the resulting preferred prominence locations would all 
be in the equatorial plane. This suggests that the polar field 
must be of mixed polarity. 

Even without the large-scale smoothing effect of a 
dipolar field component hidden in the dark polar caps, 
however, the global field topology implied by the surface 
magnetograms does not appear to be consistent with the 
large number of observed prominence locations. The coronal 
structure of both AB Dor and LQ Hya has been studied by 
extrapolating the coronal magnetic field from the surface 
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t = lo t = ti t = t, t = t3 

Figure 2. A schematic diagram of the formation of prominence-bearing loops. Initially, at t = to a current sheet is present above the 
cusp of a helmet streamer. Reconnection at in the current sheet a.t t = ti produces a closed loop ai t = t2. The stellar wind continues to 
flow until pressure balance is restored, thus increasing the density in the top of this new loop. Increased radiative losses cause the loop 
to cool and the change in internal pressure forces it to a new equilibrium at t = ts. 



magnetograms HJardine, UoUier Uameron fc Donati 2002| 
[Mclvor ct al. 20041 . While the fields close to the surface 
may be highly structured, at heights of several stellar radii 
where the AB Dor prominence systems are found, the field is 
largely smooth and resembles a highly-inclined dipole. This 
would give only two longitude bands in which prominences 
would be found. By filling these magnetic loops with hydro- 
static atmospheres, llJardine et al. 200211 have determined 
the coronal density and X-ray emission and shown them to 
be consistent in magnitude and rotational moduation with 
the observations. Comparison with simultaneous Chandra 
spectra also gives consistent results (.Hussain et al. 200511 . 
The X-ray emission typically comes from the high-latitude 
regions where the large flux concentrations are found and 
is confined within 1_R* of the surface. Above these heights, 
the pressure of the hot coronal gas is sufficient to break 
open the magnetic flux loops and escape to form the stellar 
wind. 

This presents a serious dilemna for models of promi- 
nences which require the restraining force of closed magnetic 
field lines to contain them. Such large loops cannot contain 
the ten million degree plasma of the stellar corona at the 
observed prominence heights of several stellar radii. We can 
see the same effect at work on the Sun. The distance at 
which the field becomes open can be estimated by compar- 
ing the gas pressure (assuming an isothermal, hydrostatic 
atmosphere) with the magnetic pressure of, say, a dipolar 
field. For typical solar parameters, this distance is about a 
solar radius above the surface. Above this height, the coronal 
gas is capable of pushing open the field lines and escaping 
as the solar wind. 

It seems then that the picture of compact, highly- 
structured, dense coronae that emerges from X-ray obser- 



vations and from surface magnetograms coupled with field 
extrapolations is at odds with the presence of many indi- 
vidual prominences trapped in the coronae of these stars at 
distances of 2-5 7?* above the surface. The aim of this paper 
is to present an alternative model for stellar coronae that is 
consistent with all of these observations. We suggest that the 
prominences are embedded not in the closed. X-ray emitting 
corona, but in the stellar wind. As on the Sun, large helmet 
streamers are likely to lie above active regions on the surface 
(see Fig.0. Above the cusps of these helmet streamers, cur- 
rent sheets form, separating opposite field polarities. Such 
current sheets are unstable to tearing and subsequently re- 
connection, which will form closed loops. These loops will be 
filled up by the stellar wind which will continue to fiow until 
a sufficient back-pressure develops in the loop. If the loop 
summit reaches above the co-rotation radius, then the out- 
ward pull of centrifugal forces may be enough to balance the 
inward pull of magnetic tension in the loop. If a mechanical 
equilibrium can be found, then the loop may be supported 
within the wind. 

The aim of this paper is to demonstrate that such cool 
loops, with their summits above the co-rotation radius, can 
reach an equilibrium. 



2 METHOD 

We take a simple model shown schematically in the left 
panel of Fig. |5| where the background field is closed up 
to some height ys and then open beyond that height. 
This height is traditionally referred to as the source sur- 
face HAltschuler fc Newkirk, Jr. 1969[ . We then consider the 
shape adopted by a closed loop that is embedded in this 
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Figure 3. Schematic diagram of the coordinate system used. 



background field as shown in the rightmost panel of Fig. 
|5| We consider a loop that is slender |Spruit 19 81 1 so that 
we can neglect any variation in pressure or density across 
its width. We also assume for simplicity that the loop is 
isothermal. 

At any point along this loop, pressure balance at the 
loop boundary ensures that 

where subscripts i and e refer to quantities internal and 
external to the loop respectively. The equilibrium shape of 
the loop is determined by the forces acting on it: pressure 
gradients, gravity and the Lorenz force: 



Vp- pg 



(2) 



We decompose this equation into components along (s) and 
normal to in) the field (see Fig. |^ . We choose to place our 
loops in the equatorial plane of the star where the gravita- 
tional acceleration (including the effect of rotation) is purely 
radial. In the cartesian coordinate system shown in Fig. Q 
we therefore have g = (0, g) and so we have 



9 



(3) 



{y + R* 

where y is the height above the stellar surface. Along the 
direction of the loop, force balance is very simple since the 
Lorentz force has no net component along the field. This 
determines the variation of the plasma pressure: 



Us 



P9s 



(4) 



where gs = g.B/|B| is the component of gravity (allowing 
for rotation) along the field. This reduces directly to 

dp 



dy 



= pg{y)- 



(5) 



Hence along each field line the plasma simply adjusts itself 
into a hydrostatic balance with 



Po exp 



Po exp 



m 
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(6) 
(7) 
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Figure 4. Variation with height of the plasma pressure scaled 
to its base value (top); the difference between the external and 
internal plasma pressures (middle) and the gradient of the pres- 
sure difference 9(pe — Pi)/dy = (pe — Pi)g i.e. the buoyancy (bot- 
tom) for loops that are cooler or hotter than their environment: 
Ti = lOTe (dashes), Ti = Te (dots), Ti = O.STe (dot-dashes), 
Ti = O.lTe (solid), Ti = IQ-'^Te (long dashes). For the very 
coolest loop (long dashes) the rapid drop in pressure near the 
surface and its steep rise at around can be seen as near- 

vertical lines in the top panel. The temperature of the external 
medium is Te = 2 X lO^K. Note that the pressure reaches a min- 
imum at the co-rotation radius, which here is at 1.7-R* above the 
stellar surface. 



We note that the pressure po at the base of each field line 
need not be the same, but for simplicity we choose one value 
Pos for the external field and one value poi for the internal 
field. For simplicity we take the external pressure to be in 
hydrostatic equilibrium and neglect the effect of the stellar 
wind. The balance of forces normal to the loop then de- 
termines the equilibrium shape of the loop. We define unit 
vectors along and normal to the field as: 



(1,2/') 



[l+(y')2]l/2 



(8) 
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Figure 5. The shape taken up by cool field lines at temperatures 
of 5 X lO^K (solid), 2 X lO^K (dashed) and 2 X IQ-'K (dotted). 
Within the external field, the temperature is Te = 2 X lO^K, the 
source surface = \ and we choose fc = 3 so that each helmet 
streamer of the external field spans 60° of longitude at the stellar 
surface. The pressure at the loop base is twice that of the external 
field. 



[l + (y')2]l/2' 



(9) 



where y{x) defines the path of a fiux tube and primes denote 
derivatives with respect to x. Gradients normal to the field 
can then be expressed as 



n ■ V = 



1 {> d _ d 



(10) 



If we also note that in (|5J the component of the magnetic 
tension term normal to the flux tube can be written as 



~y 



(11) 



then the component of force balance normal to the loop gives 

Bf 



dp 
By 



d ^ d\Bl 



y 



l + {y'r M 



-P9- (12) 



In taking derivatives of B^ it should be remembered that 
Using @, however, this can be 



B't = Bl + 2^i(pe 
reduced to 

(y'd _ d\ Bl 

y dx ay J 2 



Pi) 



1 + iv') 



rB 



(13) 



We choose to deflne as boundary conditions that the height 
of the loop (H) is given and also that the loop is flat at the 
top, i.e. y(0) = H and t/'(0) = 0. 



3 NATURE OF THE EQUILIBRIA 

There are two main constraints on the existence of an equi- 
librium for any set of loop parameters (poi, -Boi, T;) for a 
star of a given mass, radius and rotation rate. The first re- 
lates to the variation of the pressure along the loop. Clearly, 
the magnetic pressure inside the loop must always be pos- 
itive. From (0, however, we can see that the loop pressure 
is determined by the pressure of the external field and by 
the difference in the plasma pressure inside and outside the 
loop. While Bf. will always be positive, if (pe — Pi) is suffi- 
ciently large and negative there may be no equilibrium for 
the loop. 

We show in Fig. Q the variation of pressure with height 
for a loop whose temperature is equal to, lower than, or 
greater than its surroundings. In all cases, the pressure ini- 
tially falls (at a rate that increases as the temperature falls) 
then starts to rise at the co-rotation radius where g = and 
hence from Q 



yK_ 



/ GAL 



1/3 



1. 



(14) 



The pressure continues to rise with height until it is equal 
to its base value at a height 



1 + 



8GAL 



1 + 8 



m 



(15) 



(16) 



In Fig. we have chosen stellar parameters appropriate for 
AB Dor (Prot = 0.514 days) and so the co-rotation radius 
is at l.TRi, above the stellar surface. For clarity, the base 
pressures of the internal and external regions have been set 
equal, but changing them would only scale the curves up 
or down. For the case where the base plasma pressures are 
equal, the point j/m is also where the pressure difference 
(Pe —Pi) goes to zero and hence from Q Si = B^. The effect 
of a stellar wind would be to reduce the external plasma 
pressure and hence to increase the size of the region in which 
the buoyancy force is negative. 

We show also the pressure difference (pe — Pi) as a func- 
tion of height for these cases. The pressure difference always 
has a turning point at the co-rotation radius, but for low in- 
ternal temperatures, two other turning points are possible. 
One of these is close to the stellar surface, and one is close 
to ym- Beyond ym the pressure difference also changes sign. 

The second constraint on the existence of an equilib- 
rium comes from the balance of forces normal to the loop. 
As shown in 1131 the loop must adjust its shape (y') such 
that its Lorentz force is zero. The gradient of the magnetic 
pressure normal to the loop (the left-hand side of 11311 ) must 
exactly balance the magnetic tension (the right-hand side of 
1131 1). The equilibrium shape of the loop therefore depends 
entirely on the nature of Bf and hence, from (0 on the form 
of the external field Be and the pressure difference across the 
boundary of the loop (pe — Pi). The nature of this equilib- 
rium is most easily understood by considering the region 
close to the loop summit, y = H where the normal vector 
lies in the direction of — y. Since the magnetic tension force 
acts inwards, the gradient of the magnetic pressure must act 
outwards, i.e. we must have dBf/dy < 0. This places sig- 
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Figure 6. Variation with height of the external field strength 
squared (scaled to its base value). Values are calculated at the 
positions of the loop summits, i.e. at a; = 7r/2fc and with the 
source surface at ys = IR* and fc = 3. 



nificant restrictions on the types of loops that can reach an 
equihbrium, since pressure balance across the loop Q and 
hydrostatic equilibrium Q give 

9^ = 9^ + 2M(pe - p^)g. (17) 

Thus, at the loop summit, the downwards pull of tension is 
balanced by a combination of the pressure gradient of the 
external field and the buoyancy of the gas inside the loop. 
This buoyancy term may act outwards or inwards however 
as can be seen in the lowest panel of Fig. Q. For loops that 
are just a little cooler than their environment, the buoyancy 
changes sign only at the co-rotation radius, but for loops 
that are significantly cooler than their neighbours, two ad- 
ditional changes of sign are possible. 



4 SPECIFIC EXAMPLE 

In order to give an example of the types of loop equilibria 
that may be found, we choose a simple external field that 
is closed at low heights, then opens up beyond some given 
height 3/s. For simplicity, we choose a field that is potential 
and two-dimensional, so we ignore any field components that 
lie in the north-south direction and consider a field structure 
that lies in the equatorial plane where the effective gravita- 
tional acceleration is purely radial. We choose the form 



(18) 

-2*^= . As iUus- 



Ber. + iB,y = Bo^e2»fe- + e-2'=y= 

where z — x + iy and Bq — B^{0, 0) = 1 4- e 
trated in the left-hand panel of Fig. Q this field structure 
has a current sheet that extends from y — ys to infinity. 
Given the two constraints on the existence of a loop equilib- 
rium: 



Bf = b! + 2^i(pe 



> 



'dBf' 




dy 


y = H 



Pi)9 



(19) 
(20) 



it is helpful to consider how B^ varies with height at the 
loop summits (i.e at x = 7r/2fc). Here, B^ 



'2ky 



and so falls to zero when y 



As shown in Fig. Q, 



it rises beyond y — ys and tends to B^ — e~^^^' at large 
heights. Hence, below the source surface, the pressure gra- 
dient of the external field acts upwards and so acts to coun- 
teract the loop's own tension force. Thus loops of both pos- 
itive and negative buoyancy may be supported, provided 
the buoyancy force is not too large. Above the source sur- 
face, however, the pressure gradient of the external field acts 
downwards (in the same direction as the tension in the loop) 
although its magnitude tends to zero quite quickly. In this 
region, the only loops that can reach an equilibrium are 
those where buoyancy acts outwards and is sufficiently large 
to balance the tension force. Thus cool, dense loops are the 
most likely to find support above the source surface. As can 
be seen from Fig. ||1J, however, the range of heights over 
which buoyancy is negative at the loop summit shrinks as 
the loop temperature drops. For very cool loops it is confined 
to a region close to the value of ym- 

Clearly, for a given base pressure, the loop tempera- 
ture is crucial in determining how buoyant the loop is and 
therefore the range of loop heights that give an equilibrium. 
Another crucial factor however is the stellar rotation rate 
since this determines the height of the co-rotation radius 
where g = Q. \i the source surface is outside the co-rotation 
radius, then loop equilibria may be found out for both nega- 
tively and positively buoyant loops out to the source surface. 
Above the source surface, loops where buoyancy acts out- 
wards may also be supported. If, however, the source surface 
is inside the co-rotation radius, then there may be a range 
of heights between the source surface and the co-rotation 
radius where no solutions are possible because of the down- 
ward pressure of the external field at the loop summit. At 
the co-rotation radius itself, no solutions will be possible be- 
cause the buoyancy force is zero here. Some solutions may 
be found at larger heights however, where the pressure of 
the external field has reduced and an outward buoyancy can 
support the loop against tension. 

In order to clarify the types of loop equilibria that 
are possible, we consider a specific example. We choose the 
young solar-like star AB Dor, for which the co-rotation ra- 
dius is at 1.7R* above the stellar surface, and we choose an 
ambient coronal temperature Te = 2 x lO^K. Fig ||KJ shows 
the structure of some sample cool loops that extend above 
the source surface into the open field region. Clearly, any 
loops that form out in the open field region must have foot- 
point separations that are similar to that of the external 
arcade in which their lower sections are embedded. This im- 
mediately suggests that the range of loop geometries that is 
available for loops below the source surface is not available 
for those that form above it. 

We show in Fig. (|7|l the heights and footpoint separa- 
tions of loops at different temperatures for two cases, one on 
the left where the source surface is at a large height (Si?*) 
and where we expect to recover results similar to those of 
[Jardine fc Collier Cameron (1991}] since the opening-up of 
the field should have a minimal effect on most of the vol- 
ume of the corona. The other case, shown on the right of 
Fig. O is where the source surface is placed at (1-R*) and is 
therefore much closer to the stellar surface and affects much 
more of the coronal volume. 

If the source surface is at large heights, then most of the 
loops that may form have their summits below the source 
surface, within the closed field region. In this shown 
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Figure 7. Loop height plotted against footpoint separation for loops at temperatures of: 2.5 X lO'^K (dot-dot-dot-dash); 9.1 X lO^K 
(dot); 8.7 X lO'^K (dot-dash); 8.3 X lO^K (long dash); 7.1 X lO^K (short dash); 4.2 X 10®K (solid) . The source surface ys = 3R* (left) 
and j/s = IRi, (right). In all cases, the plasma beta at the stellar surface is 0.01, the external coronal temperature is = 2 X lO'^K, the 
co-rotation radius is at 1.7Rj, and we choose fc = 3 so that the external arcade spans about 60° longitude at the stellar surface. 



in the left-hand panels of Fig. Q , equilibrium solutions are 
available that are qualitatively similar to those found by 
[Jardine fc Collier Cameron (1991)] For a given loop temper- 
ature there are typically two types of solution: a low-lying 
case where the loop follows the path of the external field 
and a solution with a much greater summit height where 
the buoyancy of the gas within the loop significantly influ- 
ences its shape. 

If the source surface is at much lower heights, and 
particularly if it is below the co-rotation radius, then, as 
shown in the right-hand panels of Fig. Q, the nature of 
the solutions is somewhat different. Below the source sur- 
face, solutions can be found that are similar to those in 
[Jardine fc Collier Cameron (1991}] Above the source sur- 
face, only those loops whose summits have buoyancy acting 
outwards can be supported. Below the co-rotation radius, 
where g < 0, this requires that pi < p^. At the co-rotation 
radius itself, <? = and no solutions are possible. Above the 
co-rotation radius, <? > and we need to have pi > pe. In 
addition, the magnitude of the buoyancy force must be large 
enough to balance the magnetic tension at the loop summit. 
This means that only loops with a large enough value of /3 
(the ratio of plasma to magnetic pressures) will be able to 
reach an equilibrium. If /3 = then no loop solutions will be 
possible above the source surface. 

The maximum height attainable by the loop is clearly a 
function of both temperature and the position of the source 
surface. As shown in Fig. this maximum height increases 
with decreasing loop temperature. There is, however, a re- 
gion around the co-rotation radius where there are no solu- 
tions. 



5 DISCUSSION 

The aim of this paper is to demonstrate that for rapidly 
rotating stars, mechanical equilibrium exist for cool loops 
that have their summits well above the extent of the rest of 
the closed corona. The presence of co-rotating prominences 
trapped in the coronae of these stars therefore need not im- 
ply that their X-ray coronae extend out to the heights of 



beta-0.02 




Figure 8. The maximum height of individual loops as a function 
of the temperature of the loop. The external corona is at a temper- 
ature of Te = 2 X lO'^K and the star has a rotation period of 0.514 
days, giving a co-rotation radius at 1.7R*. The different symbols 
denote the height of the source surface ys which marks the extent 
of the closed corona: crosses denote loops where the source sur- 
face is above the co-rotation radius (j/s = 3R*) and stars denote 
loops where the source surface is below the co-rotation radius 

{vs = IR*)- 



2-5Rt at which prominences are observed. This new model 
for the structure of stellar coronae is therefore capable of 
explaining both the recent X-ray observations that suggest 
that these coronae are compact, with the observations of 
prominences at heights well beyond that of the X-ray emit- 
ting gas. 

The loop temperature is crucial in determining whether 
loops can be supported out in the open field region. Loops 
that extend into the open field need the outward pull of 
centrifugal forces acting on the gas inside the loop to sup- 
port the loop against its own magnetic tension. The loop 
temperature is therefore important as it determines the 
magnitude and sign of the buoyancy force. It also deter- 
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mines the plasma pressure of the gas inside the loop which 
must not be so large that the loop has to expand enor- 
mously to remain in pressure balance with its surround- 
ings. If the source surface is inside the co-rotation radius, 
then there may be a range of heights above the source sur- 
face where there are no solutions. In particular, no solu- 
tions can be found at the co-rotation radius itself since here 
the buoyancy force is zero and so cannot balance the ten- 
sion of the loop. For those cool loops where the buoyancy 
acts outwards at the loop summit, solutions may be recov- 
ered again at larger heights, just below the critical height 
ym/R* = l/2(-3 + ^/l + S[yK/R* + 1]^) where the plasma 
pressure in the loop rises to its base value. 

This critical height j/m represents the maximum height 
for all cool loops. As the stellar rotation rate increases, both 
yK and ym decrease, and the maximum loop height moves 
closer to the surface. The behaviour of the source surface 
is however much more difficult to determine, as the energy 
density of both the coronal gas and the magnetic field de- 
pend on the stellar rotation rate, but in a manner that is 
not well understood. ( |Schrijver7 DeRosa fc Title 2003} find 
that the source surface should increase rapidly with rotation 
rate, up to activity levels perhaps 10 times that of the Sun. 
The behaviour of the source surface for the very rapid rota- 
tors discussed here is not clear however. It may well be that 
observing the distribution of prominence positions is one of 
the best ways of locating the source surface. 

This simple model does not address the question of the 
stability of stellar prominences. To study the stability we 
would need to know more about their fine structure. Fu- 
ture observational studies should determine the filling factor 
of cool material in these prominences. While the equilibria 
presented in this paper demonstrate the concept that cool 
loops can be supported within the open field regions be- 
yond a star's close corona, they are not intended to model 
the structure of prominences themselves. The cool gas con- 
tained in the observed stellar prominences does not extend 
all the way to the stellar surface and the thin fiux tubes 
of this model are unlikely to be a good approximation for 
stellar prominences, unless they are composed of many fine 
cool threads. A more sophisticated model would be required 
if the full thermal and spatial structure of prominences is 
to be understood. These cool equilibria may however act as 
the starting points for prominence formation. The range of 
observed prominence positions for AB Dor (between 2 and 
5 stellar radii above the surface) is certainly consistent with 
the the values of yK = 1.7R* and ym = 4.7R*. 

If this model is correct and the prominences do indeed 
lie above the cusps of helmet streamers, then, as on the Sun, 
their positions should be correlated with the locations of the 
neutral polarity lines at the surface. The number of promi- 
nences would then be related to the number of changes of po- 
larity at the surface, which is a measure of the complexity of 
the field. On the Sun, the strength of the small-scale field rel- 
ative to the large-scale dipolar component varies through the 
solar cycle. It is smallest at cycle minimum where the field 
most closely resembles a dipole, and grows though the cycle 
as progressively more and more East- West bipoles emerge 
through the surface. It is therefore possible that if the promi- 
nence locations are related to surface polarity changes, then 
a systematic change in their number may signal changes in 
the underlying stellar field. 
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